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Classical analytic theories of the solar system indicate that it is stable, but numerical 
integrations suggest that it is chaotic. This disagreement is resolved by a new analytic 
theory. The theory shows that the chaos among the Jovian planets results from the overlap 
of the components of a mean motion resonance among Jupiter, Saturn, and Uranus, and 
provides rough estimates of the Lyapunov time (10 million years) and the dynamical lifetime 
of Uranus (10 18 years). The Jovian planets must have entered the resonance after all the gas 
and most of the planetesimals in the protoplanetary disk were removed. 
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The predictability of planetary motions was largely responsible for the acceptance of New- 
ton's theory of gravitation. In spite of this, Newton doubted the long term stability of the 
solar system. Laplace noted that the ratios \x of planetary masses M to solar mass M & are 
small (fi ~ 10~ 3 to 10~ 9 ), as are the planet's orbital eccentricities e ~ 10~ 2 and inclinations 
i ~ 10~ 2 (in radians). Neglecting terms proportional to second or higher powers of these 
quantities, Laplace showed that the motions of the planets were stable [(TJ. In this century, 
Arnold showed that for /z, e, and % of order 10 -43 , most planetary systems, in the sense of 
measure theory, are stable and undergo quasiperiodic but bounded variations in semimajor 
axis (a), e, and i for each planetj£2j. However, the values of the small quantities in our solar 



system are much larger than 10~ 43 , so the applicability of Arnold's theory is uncertain. Stud- 
ies over the last decade have been dominated by brute force numerical integration. Sussman 
and Wisdom \( 3} , and Laskar \(4\ performed numerical integrations of the planet's orbits and 
found positive Lyapunov exponents, indicating that they are chaotic. Sussman and Wisdom 
also showed that integrations of the Jovian planets Jupiter, Saturn, Uranus, and Neptune 
are by themselves chaotic. In neither case are the variations in a, e and i quasiperiodic, nor 
is it clear that they are bounded. Are the numerical results incorrect, or are the classical 
calculations simply inapplicable? 

We show analytically that the results of Laplace and Arnold do not apply to our solar 
system. The chaos seen in integrations of the outer planets arises from the overlap of the 
components of a three body mean motion resonance among Jupiter, Saturn, and Uranus, with 
a minor role played by a similar resonance among Saturn, Uranus, and Neptune. We test the 
theory using a suite of numerical integrations. The widths Aa/a of the individual resonances 
are of order 3 x 10~ 6 , so that small changes in the initial conditions of the planets can lead to 
regular motion. This explains the puzzling dependence of Lyapunov time with integration 
step size seen in the outer planet integrations of Sussman and Wisdom ^3% However, the 



uncertainties in the initial conditions, and those introduced by our numerical model, are 
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comfortably smaller than the width of the individual resonances, so our solar system is almost 
certainly chaotic. The resonance is extremely weak and hence easily disrupted. Torques 
exerted on the planets by the protoplanetary gas disk and by planetesimals were orders of 
magnitude larger than the resonant torques, so most of the planetesimals and all the gas 
must have been cleared from the outer solar system before the planets entered the resonance. 

Analytic Theory. Orbital dynamics. Planets in the solar system follow nearly Kep- 
lerian orbits*. The orbit of each planet can be thought of as consisting of three nonlinear 
oscillators, corresponding to the three spatial directions. The Kepler problem is unusual 
in that all three oscillations have the same frequency. The orbital elements were chosen to 
take advantage of this degeneracy. The angle / varies on the orbital time scale, while the 
angle w describing radial motion and the angle Q describing vertical motion are fixed. In 
the actual solar system w and Q are time dependent, with frequencies denoted by gj and Sj 
respectively. These frequencies are proportional to the mass ratios //, and are consequently 
much smaller than the mean motion n = dl/dt, the time rate of change of the mean anomaly. 
While our model contains only the Jovian planets, we label gj and Sj with j = 5, 6, 7, and 
8, corresponding to the radial order of the planets in the solar system. The mean motions n 
in units of cycles per day and the modal frequencies of the Jovian planets were determined 

°The orbits have sizes and shapes described by semimajor axis a and eccentricity e. The 
orientation of an orbit is described by the inclination i, the longitude of the ascending node 
Q, and the longitude of perihelion zu, while the location of the planet in the orbit is described 
by the mean anomaly I or equivalently the mean longitude A = I + zu. Collectively these 
variables are called orbital elements. In the Kepler problem, where a single planet orbits a 
spherical star, all the elements of the planet except the mean longitude are fixed, which is 
why the elements are useful quantities. We use the masses and a, e, and i from the JPL 
ephemeris DE200 (Table 1) 
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by numerical integration of the equations of motion (Table 2). Each planet's elements vary 
with all the frequencies s and g. For example 

ej sin wj m e 55 sm(g 5 t + f 5 ) + e 56 sm(g e t + 6e) + • • • , (1) 

where e 55 ~ 0.044, e 56 ~ 0.016, and £ 5 and £ 6 are constants. 

Resonances and chaos. A resonance occurs when two or more oscillators are coupled 
in such a way that a linear combination of their angles o = J2iPi@i undergoes a bounded 
oscillation, in which case a is said to librate. In the sum defining a, i denotes the ith 
oscillator and the piS are (possibly negative) integers. When the oscillators are not resonant, 
all possible combinations of O^s increase or decrease indefinitely, in which case a is said to 
rotate. The physical significance of a resonance is that energy is exchanged between the 
oscillators over a libration period, which is large compared to the oscillation period of any 
of the oscillators. This prolonged exchange can lead to large changes in the motion of the 
system. The orbit that divides regions of phase space where a librates from those where a 
rotates is called the separatrix. 

The other bit of dynamics needed to understand our result is the notion of resonance 
overlap. Chaos in Hamiltonian systems, of which the motions of the planets are an example, 
arises when the separatrix of one resonance is perturbed by another resonance. The extent of 
the chaos depends on the stochasticity parameter K, which is a function of the of separatrix 
width divided by the distance between resonances. If K is small, there is little chaos, but 
for K > 1 the region in the immediate vicinity of the resonances is primarily chaotic T5) 



An orbit which, at different times, both librates and rotates must cross the separatrix, and 
is therefore chaotic. Another signature of chaos is that two initially nearby chaotic orbits 
diverge exponentially with time; in our numerical work we use both diagnostics. 

Two body mean motion resonances. Two planets are said to be in a mean motion 
resonance when pid\\/dt ~ pzdX^/dt. In that case conjunctions between the planets occur 
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at nearly fixed locations in space. The designation "mean motion" is a little misleading, 
because if p\ 7^ pi there is no coupling between the (A, a) motion of two planets that does 
not involve a third degree of freedom, either the radial (w,e) or vertical (fl, i) motion of at 
least one of the planet s[[6l 



There are no two body mean motion resonances among the planets. However, there is a 
near mean motion resonance between Jupiter and Saturn; Jupiter makes five circuits around 
the sun in about the same time that Saturn orbits twice. Saturn affects the orbit of Jupiter 

through its gravity, described by the potential 

GM S 

<P=—, r> ( 2 ) 

where M$ is the mass of Saturn, rj and r$ are the position vectors of Jupiter and Saturn, 
and G is the gravitational constant. To see the resonance mathematically, we expand rj and 
in terms of the orbital elements of the two planets, keeping only the lowest order terms: 

<P = : — ^gLr( a s/ a j) ek s e J iP s ir j cos [2Aj - 5A 5 + kw s + qwj + ptt s + rfij] . (3) 

®& k,q,p,r 

The amplitudes 4>k, P , q ,r can be found in classic references \(7\ This result shows explicitly 
that the gravitational coupling between two bodies on Keplerian orbits always involves either 
(zu, e) or (Jl, i), so that at least three oscillators are affected. Symmetry considerations show 
that the integers in the argument of the cosine must sum to zero, 2 — 5 + k + q + p + r = 0, 
and that p + r must be even. To lowest order in the eccentricities and inclinations, the 
integers k, q, p, and r are non-negative and must sum to three. The strength of the coupling 
is proportional to e 3 or ei 2 , so this resonance is said to be of third order. Hence there are ten 
frequencies associated with the resonance, four involving only perihelion precession rates, 
such as 

2Aj - 5A 5 + 2wj + w s , (4) 
and six involving the precession rates of the nodal lines, including 

2Aj - 5X S + wj + tlj + tl s - (5) 
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The dot over the angles in these expressions denotes a time derivative. Each of the ten 
members of Eqn. @ is referred to as a resonant term or, sometimes, as a resonance. The 
reason for this misuse of terminology is that, while none of the frequencies associated with 
these terms in our solar system vanish, they are much smaller than the mean motions of 
Jupiter and Saturn. As a result, the resonant terms have a strong effect on the orbits of the 
two planets. 

Eighteenth century astronomers, unaware of the significance of these long period terms, 
noted a discrepancy between the predicted and observed longitude of Jupiter and Saturn. 
This discrepancy, known as the great inequality \(8)[ was finally explained by Laplace. He 
noted that the resonant terms given by Eqn. (|3]) force a periodic displacement of 21 minutes 
of arc in Jupiter's longitude and 49 minutes of arc in Saturn's, showing that the predictions 
of the law of gravitation agreed with observations of the two planets. 

The largest effect of Saturn's gravity on ej is the secular variations described by Eqn. 
(0). However, the most relevant component of Saturn's gravity for chaotic motion is described 
by Eqn. (|3|). This component forces much smaller variations in ej sinwj given by 

2,5 • ^ v-^ j.(2,5) k p-1 -q -r 

x sin[2Aj - 5A 5 + kw s + (p - + q£lj + rfig]. (6) 

The largest variation in e 2 j 5 , corresponding to k = 2, p — 1 = q = r = and 02.i,o,o ~ 9.6, has 
an amplitude of about 3.5 x 10~ 4 . Our numerical integrations yield 3.7 x 10~ 4 , consistent 
within the errors introduced by keeping only the highest order term in e. This variation in 
ej plays a central role in producing chaos among the outer planets. 

There are other two body near mean motion resonances in the solar system. Of particu- 
lar relevance here is the 7\u — A j near resonance between Jupiter and Uranus. The potential 
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experienced by Uranus is 

<P = : — ^k,'ql,r e J e v iP J ir u cos J - 7 -V + + 9^t/ + pttj + rVLjj] . (7) 

^ k,q,p,r 

(7 1) 

To lowest order (sixth) in e and z, there are 44 terms. The coefficients (Pkq pr range from 
~ 10~ 3 to ~ 10. By itself this resonance has little effect on the dynamics of the solar system. 

Three body mean motion resonances. Now consider the fact that e j sin wj varies; sub- 
stituting (H) into ([?D, we find the potential experienced by Uranus due to the non-Keplerian 
orbit of Jupiter; 

a ~ GMj _ x sr^fR \A 7 > 1 ) A 5 ' 2 ) s-p p 2 

~ ^se js ajs2^( 6 - P)06- Pl p,o,o^2,i,o,o e J e u e s 

a u p= o 

x sin[3Aj — 5A5 — 7A{7 + 7wj + p(zuu — tuj) + 2zu s ], (8) 

where ajs = aj/cis ~ 0.55 and ejs — \2 — 5(ns/nj)\ ~ 1.3 x 10~ 2 . For simplicity we have 
ignored terms involving the inclinations and kept only terms proportional to e|. This three 
body mean motion resonance is second order in the masses of the planets (both fij and us 
appear) and seventh order in e. 

Using the frequencies in Table 1, and accounting for terms involving i, we find a mixed 
e-i resonance at au ~ 19.21796 AU associated with the argument 

3Aj - 5X S - 7X V + 7g e t + 2s 7 t (9) 

We find a cluster of eccentricity resonances centered at au ~ 19.2163 AU with argument 

3Aj - 5\ s - 7\u + (2 - q)g 5 t + 7g 6 t + qg 7 t, (10) 
where < q < 2, and at au ~ 19.2193 AU with argument 

3A j - 5X S - 7Xu + (3 - q)g 5 t + 6g 6 t + qg 7 t, (11) 
where < q < 3. At the present epoch the JPL ephemeris DE200 has au ~ 19.21895 AU 
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For simplicity we have described only one type of term in the potential experienced by 
Uranus, that due to the influence of Jupiter as it moves in the potential of the Sun and 
Saturn, as reflected in ej sin zuj. There are similar contributions to the potential due to 
variations in Jupiter's other orbital elements. Furthermore, there are weaker resonances due 
to the gravity of Saturn, moving on an orbit perturbed by Jupiter, and acting on Uranus. 
Finally there are much smaller terms due to the direct perturbations of Uranus by Saturn 
and Jupiter moving on their unperturbed Keplerian orbits. 

We also find three body resonances with arguments containing 

3X S - 5\u - 7X N + 7g 6 t + (2 - q)s 7 t + qs 8 t, (12) 

at 19.2187 AU ^ a v ^ 19.2195 AU. The strength of the resonance is smaller than that of 
the resonance involving Jupiter by the ratio (HN/Hj)(tjs/ e su) ~3 x 1CT 3 . 

Overlapping resonances. The overlap of the individual resonances produces chaos among 
the outer planets. The width of a typical component resonance is 



— = § J(6 - pJ^iUXiS.Oo^-^se^e&el « 2 x 10" 6 , (13) 
a>u V 66 J s 

or Aa ~ 8 x 1CT 5 AU. We must substitute powers of either e 55 or e 56 for e 5 j~ p , depending on 
the resonant argument. This resonance width is comparable to the radius of Uranus. The 
libration period is 



To = T V ^147(6 - ^^S^o.o^&lo^/^^se^e^e^ » 10 7 years. (14) 

The precession frequencies g$ and g-j determine the distance between the component 
resonances; we find 

^ « *L ( 9 -^l) * 7 x 1(T 6 (15) 
a v 21 V n v J y ' 

The stochasticity parameter is 

Aa\ 2 



K^(n-). (16) 
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Using Eqns. (|13D and (|15D in Eqn. (|16D, we see that K ^ 1, so the motion is marginally 
chaotic. Then the Lyapunov time (the inverse of the Lyapunov exponent) is given by Ti <; T 
9} 



The chaotic nature of the system ensures that the angles in the perturbing potential (§) 
experienced by Uranus are essentially random variables. These chaotic perturbations force 
Uranus's e to undergo a random walk, exploring all values between and e cross ~ 0.5; for 
e > e cross Uranus will suffer close encounters with Saturn, and may be ejected from the solar 
system. The time for this to occur is of order \(9} 



T cross « 6 x 10 17 years (17) 

where p is the exponent of e\j in Eqn. (|8[). This estimate is uncertain by a large factor, 
possibly by one or two orders of magnitude, but it is clear that Uranus will be with us for a 
long time. The resonance closest to the actual value of a\j has p = 0. 

The discovery that the great inequality was due to the 2 : 5 near resonance between 
Jupiter and Saturn clearly had a strong affect on Laplace's views regarding determinism. 
We find it ironic that the 2 : 5 resonance plays such a strong role in producing chaos among 
the outer planets, thereby placing a limit on our ability to state the positions of the Jovian 
planets in the distant future. The fact that Laplace was the first astronomer to identify 
a three body resonance in the solar system, involving three of the Galilean satellites, only 
heightens the irony. More recently, three body resonances were shown to be responsible for 
much of the chaos seen in integrations of asteroids \(10) - 



Numerical integrations. In order to test our theory, we have integrated the equations 
of motion for the four Jovian planets using a symplectic integrator \( 12} . We chose this 



simplified model rather than including all nine planets in order to isolate the effects of the 
giant planets. To account in a crude way for the effects of the terrestrial planets, we enhanced 
the mass of the sun by the their mass, roughly a part in 6 x 10~ 6 . This ensures that the 
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location of resonances between the Jovian planets is shifted by an amount which is second 
order in this mass ratio, roughly 3 x 10~ n . This is much smaller than the uncertainty in the 
orbital elements of the planets. The orbital elements, which provide the initial conditions for 
our integrations, are known to a relative accuracy of a few parts in 10 million. For example, 
Aa/a ~ 2 x 10~ 7 (600km for Uranus) \( 1 1 J , much smaller than the size of the resonances. 



To determine whether the evolution was chaotic, we measured the Lyapunov time by 
comparing pairs of integrations in which the initial conditions differed by 1.5 millimeters in 
the x coordinate of Uranus. Using the DE200 ephemeris from JPL, we confirm the result of 
Sussman and Wisdom \(3]\ that the four Jovian planets are chaotic. We find a Lyapunov time 



of about 7 million years, consistent with our analytic result and with Sussman and Wisdom's 
result of about 5 million years, given that it is difficult to measure Lyapunov times with an 
accuracy much better than a factor of 2. 

To check the robustness of this conclusion, we have carried out integrations in which 
we varied the initial ajj in ten steps of 300km; the largest displacement was ±1, 500km, 
about twice the uncertainty in the JPL ephemeris. We employed symplectic correctors 
(12\ to ensure that the relative energy errors were less than 10 -9 , much smaller than the 



uncertainties in the initial conditions. In all these integrations we found that the orbits were 
chaotic. 

To test the prediction that the motion is marginally chaotic, we carried out various 
surveys of the dynamics of the Jovian planets in which all the initial orbital elements except 
ajj were held fixed \(14\ The integration time in each survey was 200 million years. In 
our first survey we varied the initial value of a\j in steps of 0.01AU between 18.9789 and 
19.3990AU. We found that between 19.18 and 19.399 AU more than 80% of the orbits are 
regular. Subsequently we conducted a survey in which au was varied in steps of 0.0001AU 
between 19.2141 and 19.2209. The resulting Lyapunov times are plotted as a function of 
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the initial semimajor axis au in Figure 1. We plot a point at 10 s years, corresponding to 
the integration time, if the orbit appeared to be regular. The location of our solar system 
as represented in the DE200 ephemeris is indicated by the vertical line in the figure. From 
18.9789 to 19.15 AU we find a strongly chaotic region, with Lyapunov times ranging from 
25,000 to 2 million years. Examination of the resonant argument \jj — 2X^ + vj^ reveals 
that from 18.9789 to about 19.13 AU our pseudo-Uranus is in a 1 : 2 mean motion resonance 
with Neptune. From 19.13 to 19.17 AU pseudo-Uranus is in the 7 : 1 mean motion resonance 
with Jupiter described by Eqn. (0), with a Lyapunov time ranging upwards from 100,000 
years. There are four other chaotic regions visible in Figure 1, centered at au of 19.219, 
19.26, 19.29, and 19.34 AU. All of these regions are associated with three-body resonances. 

The dynamics in the region from 19.21 to 19.225 AU (Figure 2) is controlled by the 
3Aj — 5Ag — l\y three-body resonance described in Eqn. (||). We can see the effects of 
the individual resonant terms. For ay < 19.218 AU the resonances are isolated by regular 
regions, indicating that the resonance widths are slightly smaller than the distance between 
resonances. For ajj ^ 19.218 AU nearly all the orbits have finite Lyapunov times, indicating 
that the individual resonances overlap completely. Figure 3 shows the resonant angle 3Aj — 
5\s — 7Xu + 3g$t + 6g^t (the q = case of eq. (|TT|)) f° r a u = 19.21908, about one planetary 
radius larger than the value of au used in the DE200 ephemeris. It alternates between 
libration, with a period of about 20 million years, and rotation, indicating that the orbit 
is crossing the separatrix of the resonance and confirming the chaotic nature of the orbit. 
In addition to the 3Aj — 5As — resonance, there is a resonant term involving Saturn, 
Uranus, and Neptune. Our calculations suggest that this resonance is responsible for the 
chaotic zones at 19.29 and 19.34, and plays a strong role in the chaotic zone at 19.26. 

Integrations of simpler models. In another survey, we set i = for all four Jovian 
planets and again varied au in steps of 0.01 between 18.9789 and 19.3990 AU, and in steps 
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of 0.0001 between 19.2141 and 19.2209AU. The general appearance is similar to that of 
Figure 1, showing that inclination resonances are not essential to produce chaos among the 
Jovian planets. However, the chaotic region near ay = 19.219 AU is not quite so extensive, 
and the resonances appear to be isolated, like those with ay < 19.218 AU in Figure 2. In 
yet another survey, we removed Neptune. The chaotic region at ay « 19.00 AU vanishes, 
but the chaos associated with 7Xy — Xj remains. Similarly, the chaos at ay = 19.29 AU 
and 19.34 AU is no longer present. However a chaotic region at ay = 19.219 AU and a 
very small chaotic region at 19.25 AU remain. The feature near ay = 19.219 AU is even less 
extensive than in the planar case, indicating that the effects of the 3As — 5Xy — 7Xn resonance 
are more important than the effects of inclination resonances involving Jupiter. Finally, a 
survey in which Neptune is removed and the remaining Jovian planets orbit in the same 
plane reveals no chaotic motion outside the 7Xy — Xj resonance. Apparently eccentricity 
resonances involving only the inner three Jovians do not quite overlap. They must act in 
concert either with inclination resonances or with three-body resonances involving Neptune 
to produce detectable chaotic regions. 

The epoc of resonance capture. Uranus probably did not form in the current 
resonance. Planet formation is believed to occur in disks around young stars. Evidence for 
such disks, which have lifetimes around ten million years, is now abundant, including visible, 
infrared and millimeter observations of disks around young stars |( 15 \ . The observations 



show that the disks contain both gas and particulate matter. The existence of our own 
asteroid and Kuiper belts, as well as of comets, suggest that protostellar disks contain larger 
bodies as well. Current understanding of the planet formation process suggests that planets 
migrate over substantial distances early in the history of a planetary system. Goldreich and 
Tremaine \(16} showed that torques produced by interactions between a gas disk and a planet 



can cause large scale planet migrations on timescales of tens to hundreds of thousands of 
years. Interactions between asteroids or comets and planets can also cause planet migrations 
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17 \ The recent discovery (18) of Jupiter-mass objects in short period (4 day) orbits around 



nearby stars strongly suggest that planet migration is common. 

We can compare the torques exerted on Uranus by the different processes. Jupiter and 
Saturn currently exert a resonant torque on Uranus given by 

Tm „ 100 (^oME.) „ w . (18) 

V au ) \ e JS J 

The torque exerted on proto-Uranus by the gas disk in which it formed is 

fGM Q Mjj\ 3 



T g as ~ 5.6 : [i v \i g m max (19) 

|( 16 \ In this expression the quantity m max is a measure of the gap in the gas disk produced 
by Uranus. If no such gap formed, the torque produced by the gas disk is even larger. The 
minimum mass of the solar nebula is about 10 Jupiter masses, so \i g = M gas disk/A^© ~ 0.01. 
The torque produced by interactions between Uranus and a planetesimal disk is 

Tplanetesimal ~ ( " ) ( ~T7 ) ( Ttt ) 5 (^0) 



\ djj ) \MijJ \T d 

ear > 

where is the total mass of the planetesimals that interact with Uranus, Tjj w 80 years 
is the orbital period of Uranus, and T c i ear w 10 7 years is the time for Uranus to clear the 
planetesimal disk. In units of GM & Mjj / ajj the torques are T res ~ 10 -11 , T gas ~ 10 -3 , 
and T planetesimal ~ \^ % M d jM v . The planets remain in resonance only if T res ^ T gas and 
T res ^ T p i ane tesimai- Clearly, Uranus must have been trapped in the resonance after the gas 
disk dissipated. Similarly, most of the planetesimal disk must be removed before the final 
trapping can occur. 
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Table 1. Masses, in units of the solar mass M Q , and the current semimajor axes a (in 
AU) , eccentricities e and inclinations i of the orbits of the giant planets. Data taken from 

JPL ephemeris DE200. 



Planet 


// = M/Mq 


a (AU) 


e 


i (radians) 


Jupiter 


9.548 x KT 4 


5.207 


0.04749 


0.02277 


Saturn 


2.859 x KT 4 


9.553 


0.05274 


0.04338 


Uranus 


4.355 x 10" 5 


19.219 


0.04641 


0.01348 


Neptune 


5.178 x KT 5 


30.111 


0.00820 


0.03089 
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Table 2. Orbital frequencies of the giant planets, in cycles per day. Data from our 

numerical integrations. 



Planet/mode n/2n (days L ) g (days 1 ) s (days x ) 



5 2.308 x 10~ 4 8.967 x KT 9 0.0 

6 9.294 x 10~ 5 5.965 x 10~ 8 -5.564 x 10" 8 

7 3.259 x 10~ 5 6.520 x 10~ 9 -6.328 x 10~ 9 

8 1.662 x 10~ 5 1.420 x 10~ 9 -1.460 x 10~ 9 
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Fig. 1. The Lyapunov time T L as a function of initial a v . The initial orbital elements of 
the planets are taken from DE200, except for a U: which is varied. There are chaotic two body- 
resonances at ajj ~ 19.00 and 19.12AU involving Neptune and Jupiter, respectively. There 
are also chaotic regions associated with three body mean motion resonances at ajj ~ 19.18, 
19.25, 19.29, and 19.34 AU. These involve either Jupiter, Saturn, and Uranus, or Saturn, 
Uranus and Neptune. The solid vertical line shows the actual location of Uranus. 

Fig. 2. A close up of Figure 1 around the actual value of a v . Between 19.216 and 19.218 
AU we find the individual eccentricity resonances associated with the resonant argument 
3Aj — 5As — 7A[/ + qg$ + 7g%t + (2 — q)gjt, which do not quite overlap. The resonances 
associated with the argument 3A j — 5As — 7\u + qg$ + 6g$t + (3 — q)g-jt lie between 19.218 
and 19.221 AU. 

Fig. 3. The resonant argument 3A j — 5As — 7\u + 3g 5 t + 6g e t in the case au = 19.21908, 
about one planetary radius larger than the actual value of au- The libration period is 
T ~ 20, 000, 000 years. A transition from libration to rotation occurs near 60 million 
years. A longer lasting transition from libration to rotation occurs at 160 million years. The 
Lyapunov time was measured to be about 7 million years. 



